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Abstract 

We calculate the boundary S matrix for the open antiferromagnetic spin 1/2 isotropic 
Heisenberg chain with boundary magnetic fields. Our approach, which starts from the 
model's Bethe Ansatz solution, is an extension of the Korepin-Andrei-Destri method. Our 
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1. Introduction 

In a 1 + l-dimensional theory with factorized scattering, the two-particle S matrix, 
which depends on the rapidity difference A of the two particles and which we denote here 
by -R(A), is constrained to satisfy the factorization (Yang-Baxter) equation 

i?i2(A - V) RisiX) R23{X') = R23{X') RisiX) Ri2iX - V) . (1.1) 

For particles in an n-dimensional representation of some internal symmetry group, -R(A) 
is a matrix acting in the tensor product space C" (g) C^. Moreover, R12, -R13, and R23 are 
matrices acting in C^ ^C^ ^C^, with R12 = -R (S) 1, R23 = 1 -R, etc. (See, e.g., Refs. 1, 
2 and references therein.) 

For a system with a boundary, the quantum-mechanical scattering of a particle with 
the boundary is described by a so-called boundary S matrix, which we denote here by 
K(X,^). This is a matrix acting in the space C"^, which may depend on one or more 
parameters, which we denote collectively by ^, that characterize the interaction at the 
boundary. The condition that boundary scattering be compatible with factorization is"^~^ 

Ri2{X- X') Ki{X,0 R2i{X + X') K2{X' ,0 = K2{X' ,0 Ri2{X + X') Ki{X,0 R2i{X- X') , 

(1.2) 
where 

i?2l(A)=Pl2i?12(A)Pl2, (1.3) 

and P12 is the permutation matrix in C^ (S> C"^, 

Vi2{x (^ y) = y ® X for x.yeC^. (1.4) 

Furthermore, we use the notation Ki = K ^ 1, K2 = 1® K. 

Given -R(A), Eq. (1.2) determines the general form of K{X,^) up to a scalar factor, 
which is a function of A and ^. Such boundary 5" matrices have been used to construct 



and to obtain the Bethe-Ansatz solution for open integrable quantum spin chains^"^^. For 
these apphcations, what is important is the matrix structure of K{X,^); the scalar factor 
plays no significant role. Further applications of boundary S matrices are reviewed in Ref. 
15. 

Recently, there has been interest^^"^*^ in determining the boundary 5" matrices for 
physical excitations of integrable theories. This requires determining in particular the 
scalar factor of the boundary S matrix, which contains information about possible bound- 
ary bound states, etc. 

One way to determine the scalar factor is the "bootstrap" approach: in addition to 
Eq. (1.2), one imposes on the boundary S matrix the constraints of unitarity, bound- 
ary cross-unitarity^^, and boundary bootstrap^^'^^ (which in the context of spin chains 
is known as "fusion" of boundary S matrices^). These conditions determine the scalar 
factor up to a CDD-type of ambiguity. Evidently, this procedure is a generalization of the 
bootstrap approach for determining bulk 5" matrices^. Ghoshal and Zamolodchikov^^ have 
determined in this manner the boundary S matrix for the boundary sine-Gordon model. 

In this paper, we consider an alternative approach: namely, to determine the boundary 

S matrix directly from the "microscopic" theory for the excitations - i.e., from the Bethe 

Ansatz equations. In particular, we calculate in this way the boundary 5" matrix for the 

physical excitations ("spinous") of the open antiferromagnetic spin 1/2 Heisenberg chain 

with boundary magnetic fields. The Hamiltonian of the chain is given by 

1 f^"^ 1 1 1 

n=-lJ2an-cjr,+i + —at + —(7^ \ , (1.5) 

where a are the usual Pauli matrices, and the real parameters ^± correspond to boundary 
magnetic fields. Although the bulk terms are SU{2) invariant, the boundary terms break 
this symmetry down to t/(l), which is generated by 

1 ^ 
S' = ^J2<- (1-6) 

n=l 



The Bethe Ansatz solution of this model is given in Refs. 4 and 21. 

Our approach is a generalization to the case of systems with boundaries of the Korepin- 
Andrei-Destri method 22,23^ which was devised to calculate factorized bulk S matrices for 
systems with periodic boundary conditions. A key ingredient is the quantization condition 
for a finite interval, which has recently been discussed by Fendley and Saleur^'^. Our result 
K{Xj^±) for the boundary S matrix is the diagonal matrix (see Eqs. (5. 19), (5. 30), (5.36)) 

where the scalar factor a{X,^±) is given by 

, ., . . _ r (^ + 1) r (f + 1) r (^ + i(2^^ - 1)) r (f + i(2e± + 1)) 
^ '^^^ r(f + i) r(^ + i) r(f + i(2e±-i)) r(^ + i(2e± + i)) ' ^ '"^ 

Let us compare the above result with previous field theory results. The bulk S matrix 
for the Heisenberg chain^^ coincides with the bulk S matrix for the sine-Gordon model^ 
in the limit /3^ -^ Stt, which is the 5't/(2)-invariant point. * Therefore, we expect that the 
boundary 5" matrix for the open Heisenberg chain with boundary magnetic fields (1.7), 
(1.8) should coincide with the boundary 5" matrix of Ghoshal and Zamolodchikov^^ for the 
boundary sine-Gordon model with (3'^ -^ Stt and with "fixed" boundary conditions. (For 
"fixed" boundary conditions, the field theory and hence the boundary 5" matrix are U{1) 
invariant.) We have verified that the two boundary S matrices indeed coincide, up to a 
rapidity-independent scalar factor, and with some redefinitions of variables. f 



* More precisely, the true S matrices R = VR for the two models are related by a 

unitary transformation. 

f In the notation of Ref. 18, "fixed" boundary conditions corresponds to the case /c = 0. 

Moreover, the isotropic limit is performed by rescaling ^ ^ A^ — ^ (here A = (Stt//?^) — 1), 

and taking the limit A ^ with rapidity u = —iO fixed. Our result for the boundary 5" 

matrix is recovered by making the identifications 6 = txX and ^ = 7r(^-|- — \). 
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The bootstrap result of Ghoshal and Zamolodchikov for the boundary sine-Gordon 
model with "fixed" boundary conditions has been verified using the physical Bethe Ansatz 
approach by Fendley and Saleur^^. 

We now outline the contents of the paper. In Section 2, we consider the problem of 
quantization on a finite interval. We begin by demonstrating that even the simplest of 
systems can exhibit interesting boundary S matrices. Indeed, we reanalyze an example 
which is familiar to every student of quantum mechanics - a free nonrelativistic particle - 
but with nonstandard boundary conditions. This leads to a boundary S matrix which is 
momentum dependent, and which has a pole which may correspond to a boundary bound 
state. We work out the quantization condition on a finite interval, and then discuss the 
generalization for factorized scattering. As already noted, this quantization condition is a 
key ingredient of the calculation of boundary S matrices from the Bethe Ansatz. In Section 
3 we analyze the Bethe Ansatz equations for the model (1.5) using the string hypothesis. 
In Section 4, we study the ground state. In particular, we compute the root density to 
order 1/A^, and calculate the surface energy. Although these results have previously been 
obtained ^^, the calculations serve as a useful preparation for the study of excited states. 
Section 5 is the core of the paper. There we compute the density of roots and holes to order 
1/A^ for excited states. With the help of the quantization condition, we then determine 
the boundary S matrix. We also perform a nontrivial consistency check on our result. 
There is a brief discussion in Section 6. At several points in Sections 4 and 5, we must 
approximate certain sums by integrals, being careful to keep terms of order 1/N. We derive 
an appropriate formula in the Appendix with the help of the Euler-Maclaurin formula. 



2. Quantization on a finite interval 

Nonrelativistic scattering 

As a warm-up exercise, we first consider a free l-dimensional nonrelativistic particle 
of mass m on the half-line a; > 0. Usually one demands that the wavefunction t(j{x) vanish 
at a; = 0. This is a sufficient, but by no means necessary, condition for the probability 
current j{x) = iil){xY dx'il^{x) to vanish at a: = 0. We consider instead the more general 
(mixed Dirichlet-Neumann) boundary condition 

cipix) + — ^(x) = at x = , (2.1) 

ax 

where c is a real parameter of dimension 1 /length. This boundary condition also implies the 

vanishing of the probability current at a; = 0, and is compatible with the self-adjointness 

of the Hamiltonian H = p^ /2m. Assuming energy eigenfunctions of the plane- wave form 

i)p{x) = Ae'P'' + Be-'P"" (2.2) 

(we set h = 1), we can use the boundary condition (2.1) to eliminate A in terms of B; and 
we immediately obtain 



t/jpix) = B 



p — ic 



(2.3) 



We conclude that the boundary 5" matrix is given by 



K{p) = P±^ . (2.4) 

p — tc 

The pole aX p = ic implies the existence (for c > ) of a boundary bound state with energy 

E = -cV2m. 

We remark that for c = 0, we have the Neumann boundary condition t/j'^O) = 0, and 

the wavefunction is 

i^pix) ~ cos pa: ; (2-5) 
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while for c -^ oo, we have the Dirichlet boundary condition '0(0) = 0, and the wavefunction 

is 

i^pix) ~ smpx . (2.6) 

Wavefunctions with such properties have appeared in quantum impurity problems. (See, 
e.g., Ref. 27.) 

We next consider the problem of a free particle on the finite interval —L/2 < x < L/2, 
with mixed Dirichlet-Neumann boundary conditions at both ends: 

=Fc±V(a^) + -ri^ix) =0 at a; = ±— . (2.7) 

ax 2 

Assuming again that the energy eigenfunctions are plane waves (2.2), and imposing the 
above boundary conditions, we see that p obeys the quantization condition 

e'^P^ K-{p) K+{p) = l, (2.8) 

where the boundary S matrices K^ are given by 

K^i,) = l±iS±, (2.9) 

p — ic± 

Of course, the momentum operator is not defined on a finite interval, and therefore p is 
not a momentum eigenvalue. Nevertheless, the energy is still given hy E = p^ /2m. 

Finally, if we now turn on a refiectionless potential which is localized near a; = 0, the 
quantization condition (2.8) is generalized to 

e^2p^ R{pf K-{p) K+{p) = 1 , (2.10) 

where R{p) is the S matrix for the potential. We shall now see that this formula has a 
straightforward generalization for factorized scattering. 



Factorized scattering 

Consider a system of particles with factorized scattering on an interval of finite length 
L. Each particle has some rapidity A, and in general also carries internal quantum numbers, 
such as spin. We denote the energy of a particle by £(A), and we define p{\) by the 
expression for the momentum of a particle for the corresponding system with periodic 
boundary conditions. For the case of two such particles, with corresponding rapidities Ai 
and A2, the quantization condition reads (in the notation used in Eqs.(l.l) - (1.4)) 

e'^^^^'^"- Ri2{Xi - A2) i^i(Ai , e_) i?2i(Ai + A2) i^i(Ai , e+) = 1 . (2.11) 

This quantization condition has recently been discussed in Ref. 24. As noted there, such 
a formula can be derived with the help of the Zamolodchikov-Faddeev algebra. 

3. Bethe Ansatz and string hypothesis 

The Hamiltonian 7i for the open antiferromagnetic Heisenberg chain with boundary 
magnetic fields is given by Eq. (1.5). The simultaneous eigenstates of 7Y and S^ have been 
determined by both the coordinated^ and algebraic^ Bethe Ansatz. In the latter approach, 
one constructs certain creation and destruction operators, B{X) and C(A), respectively; and 
the eigenstates are given by 

B{Xi)B{X2)---BiXM)uj+, (3.1) 

where a;"*" is the ferromagnetic vacuum state with all spins up, 

C{X)u+=0, (3.2) 

and {Aq} satisfy the Bethe Ansatz (BA) equations 

A^+i(e+-|) \ / A^+i(e--^) A fK±±y ^ yv ( Xo,-Xf3 + i \ fK±Xp+i\ 
A«-z(^+-i)y' Ua-^(^--i); [Xc-y ^}Mc. - Xp - i J [x^ + Xp-tJ ' 

« = 1,---,M. (3.3) 



The corresponding energy eigenvalues are given by 

1 ^ 1 

(plus terms that are independent of {Aq}), and the corresponding 5"^ eigenvalues are given 
by 

S' = — -M. (3.5) 

The BA equations can be written more compactly as 

M 

e2e+-l(-^a) e2^_-l(Aa) ei(Aa)^^+^ = ^ n ^2(-^a " ^p) e2(Aa + A^j) , (3.6) 

13=1 

where 

A + — 
^n(A) = -^. (3.7) 

We restrict the solutions of the BA equations as follows: 

Re(Aa)>0, Aa^O,oo. (3.8) 

Within the coordinate Bethe Ansatz approach, these restrictions can be understood by 
examining the Bethe wavefunction. Alcaraz et al?^ write the wavefunction in terms of 
momenta {/cq}, which are related to our rapidities {Aq,} by exp(z/cQ,) = (Aq + |)/(Aa — |). 
Because of the periodicity Re{kct) -^ Re{kct) + 2% of the wavefunction, one can make the 
restriction —n < Re{kct) < tt. Moreover, changing the sign of (any one) /Cq, results in 
a change of sign of the wavefunction, and so does not lead to a new independent Bethe 
state. Thus, we can make the further restriction < Re{kct) < tc. Finally, one can verify 
that the wavefunction vanishes identically for /cq, = , tt. Translating the restrictions 
< Re{koi) < TT , /cq, 7^ , TT in terms of Aq leads to Eq. (3.8) above. See also Refs. 9, 21, 
24, 26, 28. 
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We adopt the "string hypothesis" , which states that in the thermodynamic (A^ -^ oo) 
hmit, all the solutions are collections of Mn strings of length n of the form (for M„ > 0) 

A^^')-A-+i(^-j) , (3.9) 

where j = 1, ■ ■ ■ , n; a = 1, ■ ■ ■ , Mn, n = 1, ■ ■ ■ , oo; and the centers X"^ are real and non- 
negative. The total number of A variables is M = Xl^i ^-^n- 

Implementing this hypothesis in the BA equations (3.6), and then (following Takahashi^^ 
and Gaudin'^'^, and using the notation of Tsvelick and Wiegmann^^) forming the product 
n^=i o'^^^ ^he imaginary parts of the strings, we obtain a set of equations for the centers 

\n. 

d(n, 2^+-l) d(n,2C--l) 

en(A2)'^+^ n en+2C+-2z(A2) n ^n+2^-2liK) 

1=1 1=1 

oo Mm 

= i-r n n ^r.UK - A^) K^(A2 + A^) , (3.10) 

m=l/3=l 

where 



and 



-Enm(A) = e\n-m\W e|„_^|_^2(A) ■■■ e^_^^_2(A) e„+^(A), (3.11) 



,, X 1 min(n,x) iix = integer /o ion 

d{n,x) = < .„ / • . (3-12) 

' n it X 7^ integer 



In obtaining the above result, we use the lemmas 

n d(n,x) 

n^-(^^'^)= n en+x+i-2KA2) , (3.13) 

J=l 1=1 

n m 

n n ^2(A2^'^'^ - A^'-''^)) = En^{\l - A^) . (3.14) 

J=lfc=l 

Since the equations (3.10) involve only products of phases, it is useful to take the 



logarithm: 



d{n,2i+-l) d(n, 2^_-l) 

(2iV+l)Q„(A2)+ Y. ?n+2^+-2z(A2)+ Y. Qn+2^^-2l{K 

1=1 1=1 
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OO Mm 

= 27vJ: +Y.J2 [-r.UK - A^) + S_(A2 + A^)] , (3.15) 

m.=l p=l 

where qnW are odd monotonic-increasing functions defined by 

Qn(A) = 7r + zlne„(A) , -tt < q„(A) < tt , (3.16) 

and EnmW are given by 

Snm(A) = (1 — Snm)q\n-m\W + 2q'|n-m|+2 (A) H +2gn+m-2(A) + q'n+m(A) . (3.17) 

Moreover, J^ are integers in the range 

•Jmin — •Jet — '^max •> (o.ioj 

where, for the case ^± = 00, 

00 
J'^mn = 1 , Jla.=N + M^-2Y, min{m, n) M^ . (3.19) 

m=l 

There are similar formulas for the case ^± 7^ 00. (See also Refs. 21,28.) We assume that 
the integers {J^} can be regarded as quantum numbers of the model: for every set { J^ } 
in the range (3.18) (no two of which are identical), there is a unique solution {A^} (no two 
of which are identical) of Eq. (3.15). 

For definiteness, we henceforth consider the case that 2^± is not an integer, and 
therefore, d{n,2^± — 1) = n in Eq. (3.15). Moreover, for later convenience, we further 
restrict ^± > 1/2. (See Eq. (4.10) below.) 

Evidently, in order to calculate the boundary 5" matrices, we must investigate the 
low-lying excitations of the spin chain. However, as a useful warm-up exercise, we first 
consider the ground state. 
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4. Ground state 



As discussed in Ref. 21, the ground state is an 5"^ = state characterized by a fiUed 
"Fermi sea" of strings of length 1. That is, Mi = -j and M^ = for n > 1. (We consider 
the case that A^ is even.) We define the "counting" function h{X) by 

{f 
(2A + l)Qi(A) + q2u-iW + 92S-i(A) - Yl [52(A - Aj) + Q2(A + AJ 
/3=i 

(4.1 



so that the BA equations (3.15) for n = 1 are 



hiXi) = Jl, a = l,...,^. (4.2) 

We define the root density p(A) by 

In order to control terms of order 1/A, one must exercise care when passing from the sum 
in (4.1) to an integral. In particular, as shown in the Appendix (see Eq. (A. 12)), there is 
an additional term of order 1/A. Hence, the root density obeys the integral equation 

p(A) = 2ai(A) - / dX' [a2(A - A') + a2(A + A')] p{X') 
Jo 

+ — [ai(A)+a2(A) + a25+-i(A)+a2c_-i(A)] , A>0 (4.4) 

(plus terms that are higher order in 1/A), where 



an(A) 



1 dqn{X) In 1 



/oo 
rfc^ e"*'"^ e-"'"l/2 , n>0. (4.5) 

-oo 



2n dX 27r A2 + ^ 2n 

In order to solve this equation, it is convenient to extend the domain of the root density 
to negative values of A, such that it is an even function of A. That is, we introduce the 
symmetric density Ps{X) 
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This density satisfies the hnear integral equation 



1 



Ps = 2ai - a2 * ps + — (ai + a2 + 025+-! + a2^_-i) 



where * denotes the convolution 



{f*g){\)= / d\' f{\-\')g{\'). 



The solution is readily found by Fourier transforms, and is given by 



(4.7) 



(4i 



Ps{\) = 2s{\) + - [.(A) + J(A) + J+(A) + J_(A)] 



where 



.(A) 

J(A) 

J±(A) 



1 



1 



2ch7rA 27r 



duj e 



\uj\i2 



1 + e- 



(ai*s)(A) 
(a2^±-2*s)(A) 



1 

2^ 



-kl 



00 1 + e" 



1 

2^ 



(icu e' 



}\ 



-{i±-h)\^\ 



l + e- 



(4.9) 



(4.10) 



We calculate the ground state energy density e^^ = E /N with the help of Eq. (3.4). 
Again exercising care in passing from a sum to an integral, we obtain (to order 1/A^) 



jV 
2 



"N 



^E«.(^i: 



'2 



Qf=l 

00 



d\ ai(A)ps(A) - ^7^1(0) 



/ 



vac I 



(4.11) 



where e"^'^ is the energy density of the infinite chain, and the surface energy / of the chain 
is given by 

d\ ai(A) [s{\) + J+(A) + J_(A)] I . (4.12) 



TT 



/ = 2 ^ ^(0) 
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Our results for the root density (4.9) and the surface energy (4.12) agree with those of Ref. 
26. 

5. Excited states and S matrix 

As shown by Faddeev and Takhtajan^^, the isotropic antiferromagnetic Heisenberg 
chain with periodic boundary conditions has spin 1/2 excitations ("kinks" or "spinous"). 
Following their approach, we consider the case that A^ is even, and therefore the number 
of such spinous is even. Hence, the lowest-lying excited states have two spinous, and there 
are four such states: the triplet {S = 1) and singlet {S = 0) states. 

For the open spin chain, the bulk terms in the Hamiltonian are 5't/(2)-invariant; hence, 
the excitations are still spin 1/2 spinous, and there are four degenerate states with two 
spinous. However, the boundary terms break the SU{2) symmetry, and so the total spin 
is no longer a good quantum number. We now proceed to investigate the excited states, 
which we classify by their S^ eigenvalue. (Excited states of the open spin 1/2 chain have 
also been studied in Refs. 9, 28.) 

S^ = 1 state 

The 5"^ = 1 state is characterized by two holes in the sea of strings of length 1, and 
no strings of greater length. That is. Mi — -^ — 1 and M^ = for n > 1. The counting 
function h^^\X) is given by 

{f-i 
(2A + l)gi(A) + Q2C+-i(A)+Q2€_-i(A)- J2 [q2{X - XD + q2{X + X})] 
13=1 

(5.1) 
The hole rapidities Ai , A2 are given by 

/l^'^(Aa) = Ja, « = 1,2, (5.2) 
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where Ji , J2 are integers in the range 

'^min S "^1 7 "^2 S J max i \p-'^) 

and which are not in the set of integers { J^} corresponding to {A^}. 
We define the quantity a^^\X) 



,,.,,. :™ 



(5.4) 



which is the density of roots plus the density of holes. * We proceed as in the case of the 
ground state, being careful when passing from the sum to an integral (see Eq. (A. 12)), and 
we obtain the integral equation 

/•OO 

a^^\\) = 2ai(A) - / d\' [a2(A - A') + a2(A + A')] ^(^^(A') 

+ Trjail-^) + a2(A) + a2^+-i(A) + a25_-i(A) 

2 
+ J^ [a2(A- Aa) +a2(A + Aa) }, A>0 (5.5) 

a = l 

(plus terms that are higher order in 1/N). The symmetric density as (A) defined by 

^s (A)-|^(l)(_;^) ;^^0 , (5.6) 

is therefore given by 



aW(A) = 2.(A) + lrW(A), (5.7) 



where r^^'{\) is given by 



(i)(A) = s{\) + J(A) + J+(A) + J_(A) + ^ [j(A - K) + J(A + A,)] (5.8) 



* Since o^^^ depends also on the hole rapidities, the notation (t*^^''(A, Ai , A2) would be 
more accurate. However, following the usual practice, we suppress the dependence on the 
hole rapidities. 
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(plus terms that are higher order in l/N). The functions s(A) , J(A) , and J±(A) are defined 
in Eq. (4.10). 

The energy density e^ = E/N is given by 



1 ^ 



a=l 

where e]^^ is the ground state energy density (4.11), and the spinon energy e{X) is given 
by 

£(A) = 7rs(A) (5.10) 

(plus terms that are higher order in l/N). 

The spinous of the antiferromagnetic Heisenberg chain with periodic boundary condi- 
tions have the same energy, and have momentum p{X) given by^^ 

p(A) = tan"^ sh(7rA) - ^ . (5.11) 

For the open spin chain, we define p{X) by this equation. Note that 

^^2ns{X). (5.12) 

We come now to the calculation of the 5" matrix. As already noted in the Introduction, 

our approach is a generalization of the Korepin-Andrei-Destri method 22,23 pj-om Eqs. 

(5.4), (5.7), and (5.12), we obtain the relation 

£rf^_ Irfp 1 (1) 

N dX -ndX^N^^ ^^'- ^^■^'^' 

Integrating with respect to A, we obtain 

^/iW(A) = 2p(A) + ^y_ dX' e{X - X') A'\X') + n , (5.14) 

where e{x) = sign x = x/\x\. Finally, we evaluate this expression at the spinon rapidity 
Ai, recalling from Eq. (5.2) that /i(Ai) = Ji. Thus, 

ott ~ ~ IT r°^ 

— J, = 2pCXi) + -J_ dAe(Ai-A)rW(A)+7r. (5.15) 
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On the other hand, let us recaU the quantization condition (2.11) for two particles 
with factorized scattering on an interval of length L. For the spin chain, the number of 
spins A^ replaces L, and the "particles" are in fact spin 1/2 spinous of rapidities Ai and 
A2. Thus, we have 

e'^P^^'^'' Ri2{Xi-~X2) K^{Xl.^-) R21CX1+X2) Ki(Ai,^+) = 1. (5.16) 



The bulk 2-particle S matrix is 5't/(2)-invariant, and is given by 



25 



r(^ + i) r(f + 1) (AJ-zP) 
r(f + i) r(^ + i) (A 



^(A) = V.A ^^^^. .:.7. :\\ (5.17) 



where X and V are the 4x4 identity and permutation matrices, respectively. This matrix 
has the following form 



R{X) 



/a(A) 

6(A) c(A) , 

c(A) 6(A) I ' ^ ' 

V a(A), 



with a(A) = 6(A) + c(A). 

The U{1) symmetry of the Hamiltonian's boundary terms implies that the boundary 
5" matrix is of the form 

^•(^•«-("*0'%(A.«))' '^■^^' 

Our task is to explicitly determine the matrix elements a(A , ^) and (3{X , ^). 
For the S^ = 1 state, the quantization condition (5.16) implies 

e^2p(^i)^a(Ai - A2) «(Ai , ^-) a(Ai + A2) a(Ai , e+) = 1 . (5.20) 

In terms of the bulk {(j)^^') and boundary (</?) phase shifts defined by 

a(A) = e^'^'''^^), a(A,0 = e'^(^'^^ (5.21) 
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respectively, the quantization condition (5.20) becomes 

^m = 2p(Ai) + l$W, (5.22) 

where m is an integer, and the total phase shift $*^^^ is given by 

$(1) = <^(i) (Ai - As) + (/>(^) (Ai + As) + ipiXi , ^-) + V>C^i , e+) ■ (5.23) 

Comparing our Bethe Ansatz result (5.15) with the quantization condition (5.22), we 
see that the total phase shift ^^^' is given by 

/•OO 

$(^)=7r/ rfAe(Ai-A) r(^)(A). (5.24) 



It is convenient to evaluate the derivative of this expression. 



2n 



s(Ai) + J(Ai) + 2J(2Ai) + J+(Ai) + J-{Xi] 



dXi 

+ J(Ai-A2) + J(Ai + A2)] . (5.25) 

In obtaining this result, we have used Eq. (5.8) for r^^\X); and we have remembered 

to differentiate also r^^\X) in Eq. (5.24). The bulk phase shift (j)^^\X) given by Eqs. 

(5. 17), (5. 18), (5.21) satisfies 

^» = 2.J(A). (5,26) 

dX 

In view also of Eq. (5.23), we see that the boundary phase shift Lp{X ,^±) satisfies 
d(piX,C±) 



dX 



n [s{X) + J(A) + 2J(2A) + 2J±(A)] . (5.27) 



The phase shift is now readily evaluated by expressing the RHS in terms of the Fourier 
transform, and using the identity 

■rfa;^2!!L = infilii) I , Re z/ > , (5.28) 

uj 1 + e-^ V^(^) 
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and also the duplication formula for the gamma function, 



2'''-^T{z)T{z + -)=n^T{2z). (5.29) 



The result for the boundary 5" matrix element (up to a multiplicative constant) is 

"(A,^±) = e 



.^(A,e±)_r(^ + i) r(f + 1) r(^ + i(2^± -i)) r(f + i(2^± + i)) 



r(f + i) r(^ + i) r(f + i(2e±-i)) r(^ + i(2^± + i)) 

(5.30) 

S^ = -1 state 

To determine the remaining element (3{X , $,) of the boundary S matrix, we consider 
the 5"^ = — 1 state. This state is most easily described by changing the pseudovacuum. 
Hence, instead of working with the states (3.1), we work now with 

C(Ai)C(A2)---C(Am)c^", (5.31) 

where c<;~ is the ferromagnetic vacuum state with all spins down, 

i3(A)c^- = 0. (5.32) 

Sklyanin has shown^ that {Aq,} satisfy the same BA equations (3.3) as before, except for 
the replacement of ^± by — ^±. The energy eigenvalues are given by the same expression 
(3.4), and the 5"^ eigenvalues are now given by 

S^ = M-^. (5.33) 

The 5"^ = — 1 state now consists of two holes in the sea of strings of length 1, and no 
strings of greater length; i.e., M\ = ^ — 1 and M^ = for n > 1. The calculation of the 
density of roots plus holes is exactly the same as for the 5"^ = 1 state, except that we must 
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track the change ^± -^ — ^±- We find that the density is given by Eqs. (5.7) and (5.8), 
except that J±(A) is now given by 

J±(A) = -(a2^^ * .)(A) = -- / du e— ^ . (5.34) 

Moreover, the quantization condition (5.16) and the form (5. 18), (5. 19) of the 5" matrices 
imply 

e^2p(^i)A^a(Ai - A2) /3(Ai , ^_) a(Ai + A2) /3(Ai , e+) = 1 . (5.35) 

Introducing the phase shift ■i/j{X,^) by /3(A,^) = exp(z'0(A , ^)), we find that 'i/'(A,^±) 
obeys Eq. (5.27), with J± given by Eq. (5.34). We conclude that the element (3{X,^) of 
the boundary 5" matrix is given (up to a multiplicative constant) by 

/3(A , e±) = - ^^'if ~!| «(A) , (5.36) 

A-^(^±- 2) 

where a(A) is given by Eq. (5.30). This completes the derivation of the result (1.7), (1.8) 
for the boundary 5" matrix. 

S^ = states 

We have already succeeded to determine the boundary 5" matrix. Nevertheless, a good 
check on this result and on the general formalism is provided by analyzing the 5"^ = 
states. In particular, with the help of expressions for the densities, we can determine 
the "momentum" p{Xi) of a spinon in terms of its deviation from the free-particle value. 
Substituting this result, as well as the known results for -R(A) and K{X,^±), into the 
quantization condition (5.16), we obtain a consistency condition. 

To this end, we consider the 5"^ = state consisting of two holes in the sea of strings 
of length 1, and also one string of length 2; i.e.. Mi = ^ — 2, M2 = 1, and Mn = for 
n > 2. For ^± -^ 00, this is the spin-singlet {S = S^ = 0) state. The counting function 
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h^^\X) is given by (see the BA equation (3.15) with n = 1) 

h^'H^) = ^{{2N+ l)q^{X) + q2U-iW + 92C-i(A) - Yl [^^(A - A^) + q2{X + Aj 

- qi{X - Ao) - qsiX - Ao) - qi{X + Aq) - qsiX + Ao)} , (5.37) 



where Aq = Af is the position of the center of the string of length 2. The density a^'^\X) 
(defined analogously as in Eq. (5.4)) obeys the integral equation 

/•CXD 

a(°)(A) = 2ai(A) - / dX' ^X - X') + a2(A + A')] a(°)(A') 
Jo 

+ -^|ai(A) + a2(A) + a2^^-i{X) + a2^_-i(A) 

2 

+ ^ a2(A - Aa) + a2(A + Aa) 

a=l 

-ai(A- Ao) -a3(A- Ao) -ai(A + Ao) -a3(A + Ao)|, A>0 (5.38) 

(plus terms that are higher order in 1/N). We conclude that the symmetric density as (A) 
(defined analogously as in Eq. (5.6)) is given by 

af(A) = 2.(A) + lrW(A), (5.39) 

where 

r(o)(A) = r(i)(A) - MX - Ao) + ai(A + Ao)] , (5.40) 

and r^^\X) is given in Eq. (5.8). 

The position Ao of the center of the string of length 2 can be determined in terms of 
the positions Ai and A2 of the holes. We accomplish this using the BA equations (3.15) 
with n — 2: 

{2N + l)q2{Xo) + Q2C+(Ao) +?2C+-2(Ao) + ?2?_ (Ao) +?2C--2(Ao) 

2 -^ 

- 27rJf + Y, [S2i(Ao - A^) + S2i(Ao + A^)] + 2qi(Ao) + Q2(Ao) . (5.41) 

f3=l 
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(We remind the reader that, in terms of the notation in Eq. (3.15), Aq is in fact Af.) 
Passing again from the sum to an integral, and making use of the result (5.39), we are led 
to the constraint 

e2C+-2(Ao) e2^_-2(Ao) ei(Ao - Ai) ei(Ao - A2) ei(Ao + Ai) ei(Ao + A2) = 1 . (5.42) 

This is considerably more complicated than the corresponding constraint for the periodic 
chain, which implies x^^p^^^° ^^> = (^X-^ _)_ A2)/2. 
For the case ^± = 00, the constraint becomes 

ei(Ao - Ai) ei(Ao - A2) ei(Ao + Ai) ei(Ao + A2) = 1 . (5.43) 

In addition to the solution Aq = *, this constraint has the solution 



Ao = a/^ + ^ 



(Ai)2 + (A2)2 . (5.44) 



For ^± 7^ 00, Eq. (5.42) for Aq is not difficult to solve. Nevertheless, the solution is given 
by a rather cumbersome expression, which we shall not present here. 

One can verify that this 5"^ = state has the same energy (5.9), (5.10) as the S^ = 1 
state. We now express the spinon "momentum" p{Xi) in terms of its deviation from the 
free-particle value. In analogy with the 5"^ = 1 case, we have 

^Ji = 2p(Ai) + -y_ dX eCXi - X) r^'\X) + TT . (5.45) 

Substituting for A^\X) the expression (5.40), we obtain the result 

27r ~ 1 r ~ ~ 1 

2p(Ai) = — Ji + - [-$(^) + Qi(Ai - Ao) + gi(Ai + Ao)J , (5.46) 

where $*-^-' is given by Eq. (5.24). 



* The restriction (3.8) can presumably be strengthened to Re (A^) > 0, and thus, we 
discard this solution. 
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We turn now to the quantization condition. For the S^ = states, the quantization 
condition (5.16) leads to a 2 x 2 matrix equation. The two eigenvalues of this matrix 
are pure phases. (Since the matrix elements of -R(A) and K{X,^±) are known, these 
eigenvalues can be computed explicitly. However, for ^± 7^ oo, the actual expressions 
for the eigenvalues are cumbersome, and we do not give them here.) We consider the 
eigenvalue expi^^^' which for ^± ^ 00 corresponds to the spin-singlet {S = S^ = 0) 
state. t Hence, the quantization condition implies 



^i2p(Ai)Af i$(") 



(5.47) 



Substituting the result for 2p(Ai) given in Eq. (5.46) into Eq. (5.47), we obtain the 
consistency condition 

where Ai and A2 are arbitrary, and Aq is a solution of Eq. (5.42). 

This relation is readily verified for ^± = 00. Indeed, for this case, the boundary 5" 
matrix is proportional to the identity matrix (q;(A,oo) = /3(A,oo)). The 2x2 matrix 
equation which follows from the quantization condition is readily diagonalized, and the 
phase shift ^^'^' is found to be 

$(0) ^ (/,(o)(Ai - A2) + (p^^Hli + A2) + 2v?(Ai , 00) . (5.49) 

Here (j)^^\X) is the singlet phase shift 



6(A)-c(A) = e^^'°'(^) = -[^]e^^'^'(^) 



(5.50) 



where (p^^^X) is the triplet phase shift. (See Eqs. (5. 17), (5. 18), (5. 21).) Since 

$(1) = (pW(Xi - A2) + (p^^HXi + A2) + 2y:(Ai , 00) , (5.51) 



f We implicitly assume that ^^'^' is a continuous function of ^±. 
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it follows that 

e<*'°'-*'^') == e, Q(Ai - A^)) e, Q(Ai + A,)) . (5.52) 

The relation (5.48) now follows from the algebraic identity 

ei(Ai - Ao) ei(Ai + Aq) = ei Q(Ai - As)") ei Q(Ai + As)") . (5.53) 

This identity is true for arbitrary values of Ai and A2, where Aq is given by (5.44). 

We have explicitly verified the formula (5.48) also for the case ^_ = 00, ^-|- 7^ 00, and 
presumably it is true in general. This equality provides a nontrivial consistency check of 
the bulk and boundary 5" matrices and of the general formalism. 



6. Discussion 

We have demonstrated how the boundary 5" matrix for the open Heisenberg chain 
with boundary magnetic fields can be calculated directly from the Bethe Ansatz equations. 
As yet, this is the only quantum- mechanical calculation of the boundary 5" matrix. We 
have restricted our attention in this paper to the case that the bulk terms are SU{2) 
invariant and the boundary terms are U{1) invariant only for the sake of simplicity. We 
see no difficulty in extending our calculations to the case that both the bulk terms and the 
boundary terms are only U{1) invariant. However, the case where the boundary terms have 
no continuous symmetry cannot be undertaken by this approach until the corresponding 
Bethe Ansatz solution is found^^. 

The analysis of the 5"^ = states for the open chain differs significantly from that 
of the periodic chain. Indeed, for the open chain, we have seen that the position of the 
center of the string of length 2 is a complicated function of the hole positions Ai and A2 
(as well as the boundary parameters ^±); while for the periodic chain, the center of the 
string is located midway between the two holes. Naively, one might worry that this leads 
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to a breakdown of factorization. However, we have seen that factorization is maintained 
by virtue of certain nontrivial identities. We expect that a similar situation holds for the 
periodic chain with four or more holes. 

This work was supported in part by the National Science Foundation under Grants 
PHY-92 22318 and PHY-92 09978. 



A. Appendix 

We derive in this Appendix a formula for approximating a sum of the form 

. Ml 
a=l 

by an integral, to order 1/N , for large A^. Here g{X) is an arbitrary function of A which 
goes to for A ^ oo. Furthermore, {A^} are solutions of the BA equations (3.15) with 
n — 1. That is, 

h{\l) = Jl, a = l,---,Mi, (A2) 

where h{X) is the appropriate counting function, and Mi ~ N/2. There are two key points 
in our discussion: (1) we use the Euler-Maclaurin formula for approximating sums by 
integrals; and (2) we use the fact that the solution A = of the BA equations is excluded. 
The Euler-Maclaurin formula has been widely used to calculate finite-size corrections in 
Bethe Ansatz systems. (See, e.g., Refs. 26, 32.) 

Let us consider the case that there are f holes, with rapidities {Aq,} 

h{Xa) = Ja J a = l, ■■■,//. (^-3) 

We denote the union of the sets {A^} and {Aq} by {Aq}, with a = 1 , ■ ■ ■ , /, where 
/ = Ml + p. Similarly, we denote the union of the sets {J\} and { Jq} by {jo}- Evidently, 

/i(Aa)=ja=a, a = 1, ■■■,/. [AA) 
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As noted in Eq. (3.8), Aq, > 0. Nevertheless, it is useful to introduce Aq = 0. Since 
/i(0) = 0, the corresponding integer is jo = 0. 

Instead of evaluating the sum (A.l) directly, it is convenient to first evaluate the sum 

I 

N 



^E^(^-)- (^-5) 



a=0 

We approximate this sum using the Euler-Maclaurin formula (see, e.g., Ref. 33) 

V /(a + kS) = - dx fix) + - [/(a) + /(a + n6)] + ■■■. (A.6) 

Indeed, transforming to a sum over equidistant points by means of the change of variables 
h{X) = J, we obtain 



a=0 a=0 

= ^ r dj g [h -\j)) + :^[9{h -\jo)) +9{h -'(Ji))] {A.7) 
(plus terms that are higher order in 1/N). Introducing the density cr{X) 



1 dhjX) 



we obtain 



1 J2 ^(Aa) = [ 'dX a{X) g{X) + -^ [g{0) + g{Ar)] . {A.9) 

By definition of {Aq,}, 

/ Ml u 

Y, 9{K) = 9{o) + Yl ^(^a) + E aCK) . (A.io) 

a=0 a=l a=l 

Therefore, 

a=l "^0 a=l 
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Making the approximation A/ ^ oo introduces an error which is higher order in 1/A^. We 
conclude that the sum (A.l) is given by 

]^E^(^^)=/ d\a{\)g{\)--Y,9{K)-^9{Q) (A12) 

a=l "^O a=l 

(plus terms that are higher order in 1/N). The presence of the last term —^g(0) should 
be noted. 
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